Theory of Single File Diffusion in a Force Field 
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The dynamics of hard-core interacting Brownian particles in an external potential field is studied 
in one dimension. Using the Jepsen line we find a very general and simple formula relating the 
motion of the tagged center particle, with the classical, time dependent single particle reflection 1Z 
and transmission T coefficients. Our formula describes rich physical behaviors both in equilibrium 
and the approach to equilibrium of this many body problem. 
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Single file diffusion of a tagged Brownian particle, 
interacting with other Brownian particles, is a model 
for motion of a molecule or particle in a crowded one- 
dimensional environment, such as a biological pore or 
channel 0, 0] and for experimentally studied physical 
systems such as zeolites Q and confined colloid par- 
ticles [3, The confinement of the tagged particle 
by the other particles is strong and severely restricts 
the motion of the particle. The description of this sin- 
gle file motion has been of much theoretical interest 
E 1 1 1 M HO, E3, El 0- For an infinite sys- 
tem and uniform initial particle density, Harris Q and 
Levitt 0] first showed that a tagged particle exhibits 
anomalous diffusion ((xt) 2 ) ~ t 1 ! 2 , rather then normal 
diffusion ((xt) 2 ) ~ t, due to the strong many body con- 
finement effect. In recent years, two important research 
direction have begun to emerge. First, the effect of an 
external force field acting on the particles is important 
since, in many cases, pores induce entropic barriers [1] 
and are generally inhomogeneous; hence, single file mo- 
tion in a periodic potential (ll| and a confining box [l2| 
were investigated. Secondly, initial conditions have a pro- 
found effect on single file motion [ID, for example, 
particles with initial delta function distribution in space 
(rather than a uniform distribution as assumed in [a 0]) 
yield normal diffusion (l3| . This is important since if a 
potential field is acting on the particles, thermal initial 
conditions will have a Boltzmann weighting, leading to 
generally non- uniform initial conditions. 

Here we provide a general and surprisingly simple the- 
ory of single file diffusion valid in the presence (or ab- 
sence) of a potential field, V(x), as well as for thermal and 
non-thermal initial conditions. In addition, we briefly 
discuss why our main results are valid even for the case 
where the underlying dynamics is anomalous fl5| . Our 
general result reproduces those previously obtained as 
well as many new ones, by mapping the many particle 
problem onto a solvable single particle model. 

In our model, 2N + 1 identical particles with hard 
core particle-particle interactions are undergoing Brow- 
nian motion in one dimension, so particles cannot pass 
one another. An external potential field V(x) is acting 
on the particles. We tag the central particle, which has 
N other particles to its left, and N to its right. Initially 



the tagged particle is at x = 0. The motion of a single 
particle, in the absence of interactions with other parti- 
cles, is over-damped Brownian motion so that the single 
non interacting particle Green function g(x,xo,t), with 
the initial condition g(x,XQ,0) = S(x — xq), is obtained 
from the Fokker-Planck equation [la ] 
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and F(x) — —V'(x) is the force field. 

In Fig. Q] a schematic diagram of the problem is 
presented. The straight line is called the Jepsen line, 
which starts from the origin x = and follows the rule 
x(t) = vt, where v is a test velocity (7j. In the interacting 
system we label particles according to their initial posi- 
tion increasing to the right (see Fig. 1). As noticed in 
0)113 j since the Brownian particles are impermeable, ev- 
ery time a particle crosses the Jepsen line from the right 
(or left), the particle number immediately to the left of 
the line will be raised (or lowered) by one, respectively. 
Hence, the particle number immediately to the left of 
the Jepsen line defines a stochastic process decreasing 
and increasing its value +1 or —I or zero randomly. 

Following Levitt [7[ we now consider a non interact- 
ing system, equivalent to the interacting one in the large 
TV limit. We let particles pass through each other, but 
switch labels upon collision, and introduce the counter 
a(t) which increases by +1 if a particle crosses the Jepsen 
line from left, and decreases by —1 when a particle crosses 
this line from right. The event when the counter a has 
the value zero, is equivalent to finding the tagged particle 
to the left of the Jepsen line in the interacting system. 
This is the case since then the total number of crossings 
from left to right is equal the number of crossings from 
right to left. So the probability of the tagged particle be- 
ing in the vicinity of xt = vt is given by the probability 
that a = 0, i.e. by the statistics of the number of tran- 
sitions of the Jepsen line M|. In what follows we depart 
from the approach in Il8| . 

Our aim is to calculate the probability of the random 
variable a, -FV(a), and then switch vt — > xt to find the 
probability density function (PDF) of the tagged particle. 
For that we designate Pll{xq J ) as the probability that a 
non tagged particle j starting to the left of the Jepsen line 
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FIG. 1: Schematic motion of Brownian particles in a har- 
monic potential (the parabola) where particles cannot pen- 
etrate through each other. The straight line is called the 
Jepsen line, as explained in the text. The center tagged par- 
ticle is labeled 0. In an equivalent non-interacting picture, we 
allow particles to pass through each other, and we follow the 
trajectory of the particle which is at the center. 



Xq 1 < 0, is found also at time t on the left of this line. 
Plr is the probability of a particle starting to the left 
of the Jepsen line to end on the right, and similarly for 
Prr and Prl- Consider first N = 1, that is, one particle 
which starts at Xq > and a second particle which starts 
at x^ 1 < 0. Then clearly we have either a = ±1 or a — 0. 
The probabilities of these events are easily calculated, 
for example Pn=i (a = 1) = Plr(x ~ 1 )Prr(xq 1 ) is the 
probability that one particle crossed from L to R and the 
other remained in domain R. Similarly P/v=i (a = 0) = 
Pll{x q 1 )Prr{x+ 1 ) + P lr {xo 1 )P rl (x+ 1 ) and P N=1 (a = 
— 1) = Pll(x 1 )Prl(xq 1 ). Since in the non-interacting 
picture, the motion of the particles are independent, we 
can use random walk theory and Fourier analysis [19( to 
find the behavior for any N 



Pn (a) = i- r d</>Uf =1 X (</>, Wo) 



Jacf> 



(2) 



Where the structure function is 



\U,Xo j ,4) = ^Plr{x^)Pr R (4)+P L l(x^)Prr(4) 



+ Plr(x 3 )Prl(x j ) + e-^P LL (xo j )P RL (4)- (3) 

We average Eq. with respect to the initial condi- 
tions x J and Xq 3 , which are assumed to be independent 
identically distributed random variables and we find 

(P N (a)) = ±- P dcb{\(cb)) N e^ (4) 

where from Eq. ([3]) the averaged structure function is 

(AW) = 



J<f> 



'PLR){PRR) + {PLL){PRR) + {PLR){PRL)+e~^{P LL )(P RL ) 

(5) 



The averages in Eq. (|5|) are easy to find in principle, in 
terms of the Green function of the non-interacting parti- 
cle and the initial density of particles, for example 
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Il(xq) I g(x,x 0l t)dxdx 



(6) 



where the mentioned replacement vt — » xt was made. 
Here Jl(xo) is the PDF of initial positions of the particles 
which initially are at Xq < 0, similarly f R (xo) describes 
the initial conditions of the right particles. In Eq. ^ L 
is the system size, which can be taken to infinity in the 
usual way. 

Eq. ^ describes a random walk where the number of 
particles N serves as an operational time. We may use 
the Gaussian central limit theorem (CLT) to analyze this 
walk, when N — * 00. In that limit the first two moments 
of the structure function, [i\ and (12, found in the small 
4> expansion 

(A($) = l + * M i^-^20 2 + O(0 3 ) (7) 

are the only two parameters needed to determine the 
behavior of P/v(a). Defining the variance a 2 = fi2 — 
(/ii) 2 , using Eq. and the normalization condition e.g. 
(Plr) + (Pll) = 1 we find the expected result 
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o 2 = (Prr)(Prl) + (Pll)(Plr) 



(8) 
(9) 



Using the CLT we have the probability of zero crossing, 
namely a = in the N — > 00 limit 



Pn (a = 0) 



exp 



V2nNa 



(10) 



This is our first general result, valid for a large class of 
Green functions and initial conditions and thus suited for 
the investigation of a wide range of problems. 

Symmetric potential fields V(x) — V(—x), and sym- 
metric initial conditions are now investigated. The latter 
simply means that the density of the initial positions of 
the left particles, i.e. those residing initially in xo < 0, is 
the same as that of the right particles, /r(xq) — /l(—xo)- 
In this case the subscript R and L is redundant and we 
use f(xo) = fR(xo) = /t(— xq) to describe the initial con- 
ditions [20l | . From symmetry it is clear that the tagged 
particle is unbiased, namely (xt) — 0. Further, since N 
is large we may expand expressions in the exp in Eq. (I1C 
in xt, to obtain leading terms 



m = AJx T + o(x T f 



(ii) 



where we used the symmetry of the problem which im- 
plies ((Plr) - (Prl))x t =o = 0, and by definition 



AJ 



d 

dxT 



[(Plr(x t )) - (Prl (x t ))} \ Xt=q . (12) 
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Similarly 

a 2 = 2(P RR (x T )) [1-(Prr (x T ))]\ + 0(x T ). (13) 

We designate (Prr{xt))\x t =o — 7t as a reflection coef- 
ficient, since it is the probability that a particle starting 
at xq < is found at x < at time t when an average 
over all initial conditions is made 



f(x ) / g(x,x ,t)dxdx . 



(14) 



As usual the transmission coefficient T = 1 — TZ is de- 
fined through T = (Prl(xt))\x t =o- Notice that these 
reflection and transmission coefficients are time depen- 
dent single particle quantities which give useful informa- 
tion for the many body problem. Also from symmetry 
we have ^{Prl(x t ))\xt=o = --£^(Plr{x t ))\xt=o in 
Eq. (fT^|) . Hence we define j = — 8{Plr{xt)) / dxT\x T =o 
where from its definition Eq. ([6]) 



f(x o )g(0,x o ,t)dx o . 



(15) 



So j is the density of non-interacting particles at x = 
for an initial density f(xa). Using Eqs. (I10II13|) we find 
our main result the PDF of the tagged particle 



P(x T ) 



: exp 



2{(x T f 



(16) 



where ({x T ) ) = TZT /(2Nj 2 ) is the mean square dis- 
placement (MSD). 

Gaussian packet. As our first example, we consider 
particles without external forces V(x) = 0, in an infinite 
system with symmetric Gaussian initial conditions with 

a width £: /(so) = V2exp(-^/(2C 2 )))/V^ (a* > 0). 
The free particle Green function is 



g(x,x ,t) 



exp 



(x-Xq)' 

ADt 



(17) 



Using Eqs. Q14I17|1 we find the MSD of the tagged par- 
ticle [U 



N 



2D A 



-arccot 




(18) 



This solution exhibits a transition from anomalous sub- 
diffusion to normal diffusion. At short times 2Di/£ 2 <C 1 

((x T ) 2 ) ~ while at long times ({x T ) 2 } ~ For 

short times the particles do not have time to disperse, 
hence the motion of the tagged particle is slower than 
normal, increasing as i 1 / 2 since it is confined by the other 
particles in the system. 

Particles in a box. We now consider the example of 
particles in a box extending from — L to L, which was 



recently investigated using the Bethe-ansatz and numer- 
ical simulations (T^. The tagged particle initially on 
x = has N particles uniformly distributed to its left 
and similarly to its right. In this case f(xo) = 1/L for 
< xq < L. The green function g{x, Xo,t) of a single 
particle in a box with reflecting walls is found using an 
eigenvalue approach [l6| , similar to the method of solu- 
tion of the undergraduate quantum mechanical problem 
of a particle in a box. We find 



K 



E 



exp 



-D 



— — i 



4L 



(19) 



«=i,Odd 



so at t = 0, TZ = 1 since all particles initially say in 
(0, L) did not have time to move to the other side of 
box, and lim^oo TZ — 1/2 since in the long time limit 
there is equal probability for a non-interacting particle 
to occupy half of the box. The eigenvalues of the non- 
interacting particle control the exponential decay of TZ 
Eq. (|19p which in turn determines the dynamics of the 
interacting tagged particle. We also find j — 1/2L using 

Eq. (fTS")) . For short times D7r 2 i/4L < 1 we can replace 
the summation in Eq. (|19p with integration and then 
using Eq. (JT6J and p = N/L 



P(x T ) 



VP 



2 (wDt) 



1/4 



exp 



p(xtY 



WDt 



(20) 



Thus the taggedparticle undergoes a single file sub- 
diffusive process @, 0] ((x T ) 2 } ~ 2(Di) 1 / 2 /(p0F), since 
for short times the particles do not interact with walls. 
In the long time limit, the tagged particle reaches an 
equilibrium easily found using Eqs. (|16I19[) 



p(xt) 



N 

e i 



(21) 



—2 



and lim 4 _» 00 ((a;t) 2 ) = L /(2N). More generally from Eq. 
dHJ), ((x T ) 2 ) = 2TZ(1-TZ)L 2 /N which nicely matches 
and simplifies considerably, the Bethe-ansatz solution 
[12 already for N = 70 ||. 

Thermal initial conditions. If we assume that ini- 
tially the particles are in thermal equilibrium, our 
simple formulas simplify even more. If f(xa) = 
2 exp(— V(xo)/kbT)/Z where the normalizing partition 
function is 



exp 



V(x) 

hT 



dx 



(22) 



then j = l/Z. To see this, note that since we have a sym- 
metric case V(x) = V(—x) then <?(0, — xq, t) = g{Q, xq, t) 
and j = 2 f Q °° exp [-V(x )/k b T] g(0, x , t)dx /Z gives 



exp 



V(x ) 
k h T 



g(0,xo,t)dx . 



(23) 
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Where we assume that the potential is binding so a sta- 
tionary solution of the Fokker-Planck equation is reached 
i.e. free particle Eq. (fT7j) is excluded. Therefore j Eq. 
(|23p is the probability of finding non-interacting parti- 
cles on the origin, with thermal equilibrium initial con- 
ditions. Since the latter is the stationary solution of 
the Fokker-Planck operator, j is time independent and 
equal to j — exp[— V(0)/kf,T]/Z. We can always choose 
V(0) = and then j = l/Z. Using Eq. HU) we find the 
MSD of the tagged particle 



<(*t) 2 ) = 



TZTZ 2 
2N 



(24) 



Single file motion in Harmonic Potential. We 
now consider particles in an Harmonic potential field 
V(x) — mw 2 i 2 /2 where lu is the Harmonic fre- 
quency and initially the particles are in thermal equi- 
librium. The corresponding single particle green func- 
tion g(x, xq, t) describes the Ornstein-Uhlenbeck process, 
and is well known [la ] . Defining the thermal length scale 
£th = \J k}>T / muj 2 and using thermal initial conditions 
f(x ) = v^Mth exp [-(x ) 2 /2£ 2 h ], we find using Eqs. 



((^) 2 > = |f(£ 



ll 1 o 

th) < 7 jarccof 



exp (2t) - 1 



4 7T^ 

(25) 

where t — Dt/(£ t h) 2 is dimensionless time. For short 
times t <C 1 we obtain sub-diffusive behavior ((xt) 2 ) ~ 
^thV^Dt/N since then the effects of the binding field 
are negligible, while the tagged particle motion is re- 
stricted by all others leading to sub-diffusive behavior. 
For long times the tagged particle reaches an equilibrium 

((xt) 2 ) ~ 7T(£ th ) 2 /4iV. 

Equilibrium of tagged particle. In the long time limit, 
and for binding potential fields we find again a simple 
limiting behavior. First, note that lim^oo 1Z = 1/2 as 
is easily obtained from Eq. (TT4"|) and physically obvious 
for the symmetric system under investigation. Secondly, 
for any initial condition g(0,Xo,t) = l/Z when t — > oo, 
namely the Green function reached an equilibrium and 
hence we find from Eq. (fT5"|) lim^oo j = l/Z. Therefore 
in the long time limit 



Consistently, this result can be derived directly from the 
canonical ensemble, using the many body Hamiltonian 
of the system, in the large N limit, without resorting 
to dynamics as we did here. The long time limit of our 
examples of particles in a box Eq. (|2ip and in a harmonic 
field Eq. (|25|) are easily recovered from the more general 
Eq. ([261). 



Z 2 

lim ((xt) ) ~ • 



(26) 



In the examples analyzed so far, we considered sym- 
metric external fields, and initial conditions, which yield 
no drift (xt) = 0. Our general results Eqs. (|10|16|) can 
be extended to cases where the drift is non-zero, then ex- 
pansions Eqs. (|11I12I13P include additional terms which 
make the final expression slightly more cumbersome and 
hence suited for a longer publication. Another interest- 
ing extension of this work is the case where the particle 
dynamics, in the absence of interactions with each other, 
is anomalous as was considered in [3, for the case 
of sub-diffusing particles [l5[ ■ These results will be ad- 
dressed in the future. 



To conclude, we have mapped the many body problem 
of interacting hard core Brownian particles to a single 
particle problem where calculation of the reflection coef- 
ficient 1Z and j yield the motion of the tagged particle. 
In the large N limit, information on the motion of the 
tagged particle is contained in the single particle Green 
function g(x,xo,t) which can be calculated with known 
methods. Rich physical behaviors emerge, which depend 
on the initial distribution of the particles and the force 
field. For particles in a confining field, a transition from 
sub-diffusion behavior to an equilibrium value is found 
to be controlled by the eigen modes of the single par- 
ticle green function, thus the approach to equilibrium 
is fast in this sense, even though the diffusion is slower 
than normal. In contrast, for Gaussian initial conditions 
a transition between sub-diffusive to normal diffusion is 
found, which is the expected transition for a wide vari- 
ety of well behaved initial conditions in the absence of 
a force. The most rewarding part of this work are the 
simple equations we find for the single file problem. 
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